如何看懂vcf格式文件 | 全外专题

1. VCF格式说明

VCF格式：Variant Call Format，用于记录variants (SNP/INDEL)的文件格式，关于其说明，详见：

http://gatkforums.broadinstitute.org/discussion/1268/how-should-i-interpret-vcf-files-produced-by-the-gatk 。

VCF是用于描述SNP，InDel和SV结果的文本文件。VCF格式在GATK软件中得到很好的支持。

2. VCF文件结构

VCF文件分为两部分内容：以“#”开头的注释部分和没有“#”开头的主体部分。（注释部分和主题部分）；注释部分有很多对VCF的介绍信息；主体部分包含10列数据。主题部分每一行代表一个variant的信息。 

主体部分10列的范例：

chrM 150 . T C 7766.77 PASS

AC=2;AF=1.00;AN=2;DP=199;ExcessHet=3.0103;FS=0.000;MLEAC=2;MLEAF=1.00;MQ=49.78;QD=32.91;SOR=0.904

GT:AD:DP:GQ:PL 1/1:0,175:175:99:7795,531,0

3. 主体部分的10列分别代表的意义

4. VCF文件的基因型信息

VCF文件的主体部分的第9列是基因型信息的多个标签，这些标签之间以冒号分割，其对应的值位于第10列，同样以冒号分割，表示第一个样品的基因型结果。若有多个样品，则VCF文件超过10列，且第10列后的每一列表示一个样品的基因型结果。第9列各个标签的意义展示如下：

GT： Genotype

样品的基因型（genotype），两个数字中间用‘/’分开，这两个数字表示双倍体的sample的基因型。0表示样品中有ref的allele（可初步理解为和ref的碱基相同，即和REF相同）；1表示样品中的variant的allele（可以理解为和variant变异后的碱基相同，即和ALT相同）；2表示有第二个variant的allele（和ALT的第二种碱基相同）对于SNP是指单个碱基类型相同而对于Indel是指碱基类型及个数均相同。

因此根据GT的结果得出以下结论：

0/0表示sample中该位点为纯合位点，和REF的碱基类型一致

0/1表示sample中该位点为杂合突变，有REF和ALT两个基因型（部分碱基和REF碱基类型一致，部分碱基和ALT碱基类型一致）

1/1表示sample中该位点为纯合突变，总体突变类型和ALT碱基类型一致

1/2表示sample中该位点为杂合突变，有ALT1和ALT2两个基因型（部分和ALT1碱基类型一致，部分和ALT2碱基类型一致） 

AD和DP

AD(Allele Depth)为sample中每一种allele（等位碱基）的reads覆盖度，在diploid（二倍体，或可指代多倍型）中则是用逗号分隔的两个值，前者对应REF基因，后者对应ALT基因型；

DP(Depth)为sample中该位点的覆盖度，是所支持的两个AD值（逗号前和逗号后）的加和，例如：

1/1:0,175:175—GT:AD(REF),AD(ALT):DP 

0/1:79,96:175 

1/2:0,20,56:76 

这里的三种类型对应的DP值均是其对应的AD值的加和，1/1的175是0+175，0/1的175是79+96，1/2的76是0+20+56。

GQ（基因型存在的概率）

基因型的质量值（Genotype Quality）。Phred格式（Phred_scaled）的质量值，表示在该位点该基因型存在的可能性；该值越高，则Genotype的可能性越大；计算方法：Phred值=-10*log(1-P)，P为基因型存在的概率。（一般在final.snp.vcf文件中，该值为99，为99时，其可能性最大）。

PL（likelihood genotypes）

指定的三种基因型的质量值（provieds the likelihoods of the given genotypes）；这三种指定的基因型为（0/0，0/1，1/1），这三种基因型的概率总和为1。该值越大，表明为该种基因型的可能性越小。Phred值=-10*log(P)，P为基因型存在的概率。最有可能的genotype的值为0。

例如：

0/1:889,0,216 

0/1:94,0,940 

1/1:269,18,0 

1/1:580,54,0 

1/2:3365,1522,1357,1842,0,1706 

1/2:307,190,178,117,0,104 

（0/0型3个数字，第一个为0

0/1型3个数字，中间为0

1/1型3个数字，最后一个为0

1/2型6个数字，倒数第二个为0）

5. VCF文件第8列信息

第8列的信息包括18种，都是以“TAG=Value”，并使用分号分隔的形式，其中很多的注释信息在VCF文件的头部注释中给出，下面对常用的TAG进行解释：

AC，AF和AN

AC（Allele Count）表示基因型为与variant一致的Allele（等位碱基）的数目；AF（Allele Frequency）表示Allele的频率，AF值=AC值/AN值；AN（Allele Number）表示Allele的总数目。比如：对2个sample的双倍体进行测序，则AN值为4。若REF上位点碱基为A，而2个sample在该位点分别为A/T和T/G，则AC值为2，1；AF值为0.50，0.25。AC：variant数目，AF：频率，AN：总数目

DP（reads覆盖度）

表示reads被过滤后的覆盖度。

FS

FisherStrand的缩写，表示使用Fisher’s精确检验来检测strand bias而得到的Fhred格式的p值，该值越小越好；如果该值较大，表示strand bias（正负链偏移）越严重，即所检测到的variants位点上，reads比对到正负义链上的比例不均衡。一般进行filter的时候，推荐保留FS<10~20的variants位点。GATK可设定FS参数。

ReadPosRandSum

Z-score from Wilcoxon rank sum test of Alt vs. Ref read position bias.当variants出现在reads尾部的时候，其结果可能不准确。该值用于衡量alternative allele（变异的等位基因）相比于reference allele（参考基因组等位基因），其variant位点是否匹配到reads更靠中部的位置。因此只有基因型是杂合且有一个allele和参考基因组一致的时候，才能计算该值。若该值为正值，表明和alternative allele相当于reference allele，落来reads更靠中部的位置；若该值是负值，则表示alternative allele相比于reference allele落在reads更靠尾部的位置。

进行filter的之后，推荐保留ReadPosRankSum>-1.65~-3.0的variant位点。

MQRankSum

该值用于衡量alternative allele上reads的mapping quality与reference allele上reads的mapping quality的差异。若该值是负数值，则表明alternative allele比reference allele的reads mapping quality差。进行filter的时候，推荐保留MQRankSum>-1.65~-3.0的variant位点。

相关阅读

作者：genome_denovo（https://blog.csdn.net/genome_denovo ）

原文来源:  https://blog.csdn.net/genome_denovo/article/details/78697679

侵删